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THE  HELICOPTER  PILOTAGE  MODELLING 
THROUGH  THE  NUMERICAL  METHOD 


Kazimierz  Szumanski 

This  work  deals  with  the  numerical  simulation  of  the  helicopter 
maneuvers  which  are  included  in  the  flight  assignment.  It  also  gives 
tne  principal  foundations,  computational  model  as  well  as  an  example 
of  the  analysis  of  the  helipopter  interrupted  take-off. 

Besides  that  the  work  also  includes  the  suggestion  of  the 
employment  of  the  above-mentioned  method  in  other  areas  which  are 
connected  with  the  evaluation  of  the  dynamics  of  similar  objects  be- 
ing operated  by  the  man. 

1.  Introduction 

The  purpose  of  this  study  is  to  find  through  the  computational 
method  a solution  of  the  course  of  the  flight  assignment  fulfillment 
which  is  composed  of  the  helicopter  maneuvers'  set.  Needs  for  this 
type  of  evaluation  come  from  the  necessity  to  deepen  the  analysis 
of  physical  phenomena  taking  place  on  the  flight  different  stages. 
Especially  it  is  justified  for  helicopters  due  to  their  multiparame- 
tric  steering  process  which  is  difficult  for  a simultaneous  evalua- 
tion at  a given  moment.  An  attempt  has  been  made  to  create  a mathe- 
matical projection  of  the  helicopter  pilotage  technique  in  the  field 
of  the  situation  evaluation  and  this  of  the  manner  of  the  maneuver 
performance . 


The  method  which  has  been  shown  on  the  example  of  the  evaluation 
of  the  helicopter  flight  dynamics  may  be  also  employed  for  solving 
the  problems  connected  with  the  maneuver  set  of  a plane  or  glider 
(e.  g.  a number  of  acrobatic  tricks)  as  well  as  for  the  gliding 
flight  modelling  (e.  g.  the  "dolphin"  technique)  or  for  the  evaluation 
of  other  steerable  machines  belonging  to  the  similar  structure,  that 
is  the  "man  operator-machine  structure  (e.  g.  car  driving).  The 
helicopter  motion  modelling  is  becoming  justified  because  of  an 
urgent  necessity  of  projecting  and  specifying  the  actual  course  of 
phenomena  arising  during  the  flight  of  the  helicopter  which  in  general 
is  an  object  purposely  steerable.  In  this  case  there  is  a necessity 
in  the  values  which  are  directly  being  solved  in  the  problem,  (e.  g. 
the  motion  parameters)  as  well  as  in  the  secondary  information  which 
can  be  found  during  the  problem  solving,  e.  g.  the  helicopter  elements' 
loading  as  a function  of  time  for  an  assumed  flight  profile. 

One  can  list  the  following  fields  in  which  the  need  for  the 
above-mentioned  type  of  information  becomes  essential: 

1.  the  resistance  analysis  (connected  both  with  the  ultimate  resis- 
tance and  with  the  fatigue  one)  for  propulsion  (e.  g.  transmissions; 
the  loading  range  is  necessary)  and  for  the  fuselage  elements; 

2.  the  reliability  analysis  - the  load  distribution  as  a function 

of  time  for  the  typical,  standard,  medium  and  extreme  flight  profiles; 

3.  designing  the  helicopter  equipment  - looking  for  the  response 
to  the  operational  process  as,  for  instance  in  the  fuel  system  (the 
fuel  consumption  which  depends  on  the  engine  power  as  well  as  on 
the  equipment  loading)  as  a function  of  time; 

4.  testing  during  the  flight  becomes  essential  for  optimization  of 
maneuvers,  for  analyzing  the  extreme  case,  for  evaluation  of  the 
maneuvers  which  are  not  typical  or  dangerous  during  their  evaluation 
through  the  test  in  the  flight; 
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5.  For  designing  new  helicopter  constructions  - for  gaining  the 
better  analysis  of  the  helicopter  flight  characteristics  on  first 
stages  of  the  project. 

The  computing  method  is  based  on  the  digital  simulation  of 
the  helicopter  dynamics.  It  is  justified  by  the  problem  complica- 
tion, nonlinearity  and  a complicated  form  of  the  interdependences, 
connections  between  the  degrees  of  freedom  of  the  helicopter  in  space 
as  well  as  the  necessity  maintaining  high  flexibility  of  the  model 
because  of  various  character  of  the  phenomena  which  are  being 
modelleu . 

A process  which  is  fixed  in  this  way  enables  to  realize  the 
following  dynamic  tasks: 

- it  serves  for  the  estimation  of  the  purposeful  maneuver  (the  closed 
system);  here  the  system  with  the  feedback  comes  into  being.  This 
system  is  supposed  to  project  three  types  of  maneuvers:  a)  the 
extreme  maneuver  (rough.  Incorrect);  b)  the  typical  maneuver;  c) 

the  optimal  maneuver; 

- it  can  serve  as  a simulator  (the  open  system)  for  the  pre-set  course 
of  the  steering  function;  it  has  to  respond  as  a course  of  changes 

in  the  parameter  set; 

- as  the  a)  problem  restriction  there  is  a possibility  to  project  the 
flight  with  the  autopilot  whose  effectiveness  is  limited  for  the 
pre-set  nominal  trajectory; 

- it  has  to  embrace  also  the  classical  methods  of  the  dynamics  evalua- 
tion such  as  uncontrolled  whirling  motions  during  the  linearization  of 
the  motion  parameter's  deviations,  e.  g.  for  the  evaluation  of  the 
free  motion  frequency  or  shape. 
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In  works  No.  [1],  [2],  and  [3]  which  are  mentioned  in  the 
references,  the  steering  by  means  of  devices  is  based  on  the  princi- 
ple of  autopilot,  whose  general  task  is  compensating  the  deviations 
of  the  actual  parameter  set  from  the  pre-set  one.  The  deviation 
evaluation  takes  place  at  the  actual  moment,  and  then  the  set  of  the 
nominal  parameters  must  be  determined  at  any  moment.  In  works  [2] 
and  [3]  it  is  the  nominal  landing  trajectory  and  its  automatic  steer- 
ing is  based  on  the  achievement  of  proper  landing  by  means  of  con- 
tinuous correction  of  the  flight  trajectory  according  to  the  pattern. 

The  possibility  of  deviation  compensation  depends  on  so-called  auto- 
pilot effectiveness  and  the  complexity  of  the  pre-set  flight  para- 
meters depends  on  the  autopilot's  perfection.  As  an  example  of  the 
pre-set  (pattern)  parameters  one  can  mention  simple  parameters,  such 
as:  maintaining  the  course,  flight  height,  descent  trajectory,  etc. 
or  complex  parameters,  such  as:  the  turning  performance,  or  automatic 
landing  [2]  and  [3]  where  a perfect  autopilot  is  indispensable.  If 
in  the  computer  science  we  can  distinguish  such  terms  as  the  zero  j 

condition  and  deviation  evaluation,  then  in  this  work  the  problem 
lies  rather  in  the  field  of  looking  for  the  zero  conditions.  When 
the  autopilot  is  being  employed  the  problem  falls  in  the  category  of 
deviation  problems.  The  difference  between  the  flight  simulator  and 
the  system  which  has  been  proposed  in  this  work  lies  in  the  simulator 
itself,  because  it  is  an  open  system  responding  to  the  pre-set  steer- 
ing impulses.  Therefore  it  belongs  to  the  component  parts  of  the 
proposed  model. 

The  benefits  of  the  employment  of  the  proposed  model  in  compari- 
son with  the  methods  emphasized  in  the  literature  on  this  subject 
are  following:  1)  the  possibility  to  solve  completely  the  flight  ex- 
treme and  non-typical  forms.  This  goal  in  the  helicopter  problem 
is  especially  essential,  since  fragmentary  simplified  analyses  pre- 
sently applied  are  not  sufficient;  2)  taking  into  consideration 
the  unification  of  the  model  and  method  and  improvement  and  precision 
of  solution  of  the  typical  maneuvers  and  auxllliary  problems  of 
dynamics  of  the  helicopter  flight. 


2.  The  Model  of  the  Helicopter  Dynamics 
in  the  Unsteady  Flight 

Any  flight  assignment  can  be  presented  as  a flight  profile 
consisting  of  the  set  of  the  N maneuvers  (the  examplary  flight  pro- 
file is  on  Figure  1).  Every  nth  maneuver  can  be  put  apart  as  a 
typical  flight  phase,  during  which  from  the  t moment  to  the  t 
moment  the  whole  principle  is  in  operation  or  the  right  to  perform 
this  maneuver  is  acting.  During  the  acceleration  it  will  be,  for 
instance,  maintaining  a permanent  height  and  the  fuselage  gradient 
(which  does  not  exceed  10-15°);  for  the  turn  it  will  represent 
maintaining  the  fixed  flight  beam  and  flight  height;  for  the  climb 
it  will  represent  the  accomplishment  of  the  pre-set  height  and  climb 
horizontal  speed  when  the  complete  power  excess  is  being  used.  All 
maneuvers  begin  at  the  n point,  where  all  parameters  of  the  helicopter 
motion  are  defined.  These  parameters  are  identical  to  the  final 
parameters  of  the  preceding  maneuver  in  the  chain  of  the  flight  phases 
which  are  defined  by  the  flight  assignment  (Figure  1).  In  the  given 
flight  assignment  or  during  the  process  of  solving  the  isolated 
n-maneuver  there  must  be  given  certain  typical  parameters  of  the  end 
of  the  maneuver  at  the  n+1  point;  these  parameters  are  being  consid- 
ered as  the  goal  parameters  of  the  pre-set  flight  phase.  In  the  phase 
of  giving  an  impetus  it  will  be  the  speed  arising  at  the  end  of  giving 
the  acceleration,  e.  g.  Vxl=40  per  hour,  for  the  climb  the  typical 
characteristic  will  be  represented  by  the  height,  e.  g.  Z2=10  m and 
the  horizontal  speed  Vx2  is  equal  to  80  km  per  hour,  etc.  The 
remaining  parameters  of  the  helicopter  motion  will  be  computed  on 
results  from  the  simulation  model  of  the  flight  dynamics  in  the  process 
of  solving  the  chosen  n-maneuver. 

The  n-maneuver  solving  is  based  on  the  definition  of  changes’ 
course  as  a function  of  time  of  the  linear  and  angular  parameters 
of  the  helicopter  motion  (location,  speed,  acceleration)  as  well  as 
on  the  steering  process  being  realized  by  means  of  the  wobble  plate 
of  the  main  rotor  (the  total  and  intermittent  pitch)  and  of  the  tail 
propeller  (the  intermittent  pitch). 


Figure  1.  The  exemplary  flight  profile 
consisting  of  the  N-maneuver  set  which 

correspond  to  individual  flight  phases.  j 

Key:  (l)  Accelerating;  (2)  ascending;  ’ 

(3)  turn;  (4)  manuever  n;  (5)  flight  in 

turbulen  atmosphere;  (6)  descending;  j 

(7)  landing.  , 
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The  special  solving  method  has  been  accepted.  Its  name  is  the 
"step  by  step"  solution  method  where  a short  period  of  time  - At  was 
accepted  (e.  g.  0.25  second)  which  secures  an  additional  division 
of  the  maneuver  on  the  I amount  of  periods.  The  "distance"  between 
the  above-mentioned  periods  of  time  is  At. 

During  the  numerical  simulation  of  the  chosen  maneuver  (the  simu- 
lation concerns  the  way  of  solving  the  flight  dynamics  as  well  as 
the  helicopter  pilotage  simulation)  it  is  prohibited  to  surpass 
certain  limits,  into  which  because  they  are  common  for  most  maneuvers 
one  can  include: 

- available  engine  power; 

- design  limits  imposed  on  the  wobble  plate  in  relation  to  shaft; 

- limitation  of  the  rotor  carrying  capacity  because  of  jet  separa- 
tion; 


Additional  limitations  arising  due  to  the  specifics  of  the  man- 
euver performed  can  be  following: 

- limitations  of  the  helicopter  gradient  and  these  of  acceleration 
(.the  requirement  of  the  flight  comfort); 

- the  limitation  of  steering  rate  by  means  of  the  steering  system 

(e.  g.  steering  dampers  which  are  used  sometimes  for  the  tail  propel- 
ler or  the  tempo  limitation  which  is  expressed  numerically  for  the 
evaluation  of  the  "gentle-abrupt"  steering  way); 

The  maneuver  can  be  realized  in  two  ways: 

1)  when  the  principle  of  the  flight  during  the  maneuver  does  not 
change  and  the  end  maneuver  parameters  can  be  obtained  only  through 
continuous  correction  of  the  flight  parameters.  As  an  example  in 
the  case  of  acceleration  the  height  and  gradient  are  kept  stable 


f 

I 


and  acceleration  ends  up  at  the  moment  of  acquiring  the  pre-set  speed. 
The  same  thing  happens  during  the  climb  to  the  pre-set  height  with 
acceleration  lasting  to  the  moment  of  the  established  speed  being 
achieved.  Continuous  correction  of  the  a and  a accelerations' 

Z X 

impulses  is  being  carried  out  in  such  a way  that  an  available  excess 
power  could  be  used  up  and  this  correction  makes  it  possible  to 
obtain  the  pre-set  goal  parameters  when  the  remaining  parameters 
change  only  according  to  the  equations  of  motion  of  the  helicopter. 


2)  when  the  flight  principle  changes  in  the  isolated  n-phase  of  the 
flight  and  there  is  no  possibility  oi  dividing  the  maneuver  into 
two  parts  because  of  their  linkage.  The  moment  of  the  2nd  part's 
coming  out  depends  on  the  flight  actual  parameters  according  to  the 
principle  of  the  flight  realization  in  the  first  part.  The  above- 
mentioned  way  can  be  applied  to  solving  all  maneuvers  which  require, 
for  instance  initial  acceleration,  and  later  braking  the  motion,  for 
few  parameters  which  are  strictly  determined  as  goal  parameters.  As 
an  example  one  can  mention  the  vertical  climb  over  the  pre-set  height,  j 

from  hovering  at  the  z^  height  to  hovering  at  z^  height.  The  heli- 
copter acceleration  takes  place  before  the  climb  begins  being  braked. 

Another  example  is  the  autorotation  landing.  The  helicopter  flies 
to  a certain,  difficult  to  determine  in  advance,  height  by  gliding 
flight,  then  the  process  of  landing  takes  place  through  the  kinetic 
energy  of  the  progressive  motion  as  well  as  through  the  inertial 
energy  of  the  main  rotor  system  in  order  to  reduce  the  descent.  Then 
during  the  maneuver's  solving  besides  continuous  "directing"  towards 
the  goal  at  the  moment  of  flight  realization  according  to  the  1st 
principle,  it  is  necessary  to  carry  out  a continuous  computation  of 
"the  braking  distance"  according  to  the  2nd  principle  in  order  to 
obtain  in  this  case  the  established  goal  parameters. 

The  steering  simulation,  enabling  to  select  (at  the  ti  time)  the 
temporary  flight  parameters  on  the  grounds  of  the  set-up  goal  para- 
meters, realizing  a maneuver  through  the  direction  effect  as  well  as 
continuous  correction  of  the  flight  with  taking  into  account  the 
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actual  and  expected  possibilities  and  limitations  of  the  helicopter  -j 

system,  is  called  the  motion  forseeing  system.  This  system  simulates  ] 

the  pilot's  actions  and  has  to  be  determined  mathematically  in  details 

if  it  must  be  used  for  the  purpose  of  computations.  In  reality  the 

pilot  realizes  this  system  in  an  intuitive  way  depending  on  the  kind  ( 

of  training,  maintaining  in  approximation  the  maneuver  principal 

parameters.  In  the  computing  program  (the  scheme  is  on  Figure  2)  the 

goal  and  principle  of  the  realization  of  the  n-maneuver  are  coded  in 

the  main  block  which  controls  the  PR^.  The  remaining  computations  | 

which  are  common  for  all  maneuvers  are  coded  in  proper  procedures. 

These  procedures  are  following: 

1 

- computation  of  the  necessary  power  designated  as  PN; 

- computation  of  the  helicopter's  temporary  equilibrium  in  the  space 

RS,  1 

- computation  of  the  N rotor's  critical  carrying  capacity;  \ 

- computation  of  the  TET  wobble  plates'  location,  ' , 

- evaluation  of  the  parameter  change  from  the  i-point  to  the  i+1 
point  in  the  phase  of  one  KR  step, 

- control  over  the  achievement  of  the  T goal, 

- computing  the  steering  derivatives  in  order  to  evaluate  the  D 
steering  reserves. 


- evaluation  of  the  limit  time  or  limit  location  for  the  previously 
specified  second  type  F maneuver  (very  often  It  is  included  in  the 
PRn  program) . 

The  next  part  of  the  section  contains  a detailed  discussion  of 
blocks  for  the  computational  procedure  given  in  the  schematic 
diagram  (Figure  2). 


Figure  2.  The  scheme  of  the  numerical  model  organization 
for  the  dynamics  simulation  of  the  helicopter  flight. 

Key:  (1)  Calculation  of  power  PN;  (2)  data  input;  (3) 

calculation  of  lift  of  rotor  N;  (4)  Calculation  of  equili- 
brium RS;  (5)  calculation  of  control  position  TET;  (6)  pro- 
cedures; (7)  calculation  of  control  derivatives  D;  (8) 
maneuver  n;  (9)  maneuver  n+1;  (10)  estimation  of  change  in 
parameters  of  motion  during  step  KR;  (11)  maneuver  N;  (12) 
monitor  accomplishment  of  task;  (13)  yes;  (14)  no;  (15) 
results  of  motion  parameters  in  step  i;  (16)  end  of  calcu- 
lation . 
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2.1.  Characterization  of  individual  Blocks  of 
the  Maneuver  Computational  Simulation 

2.1.1.  The  WE  (or  input)  block  or  the  input  datum  block 

It  is  necessary  to  submit  the  estimated  principal  parameters 
of  the  beginning  and  end  of  the  pre-set  typical  segment  of  the  man- 
euver as  well  as  the  principal  requirements  (restrictions)  or  also.' 
how  to  realize  the  transition  (the  motion  parameters)  between  the 
initial  and  final  point  of  the  set-up  segment.  More  often  it  will 
be  the  symbol  referring  to  the  proper  procedures. 

Therefore,  dealing  with  the  input  data  one  should  take  into 
consideration  the  helicopter  location  in  space: 

- the  Xy,  y^,  Zq  "co-ordinates  of  the  center  of  the  helicopter  mass 
and  the  angles  (angular  displacement)  <p  , (roll,  .pitch,  and  yaw); 
component  linear  speeds  V^,  V , which  were  taken  from  the  data 
concerning  the  model;  limiting  values  for  the  rate  of  movements  with 
the  controls,  reduced  to  the  wobble  plate  i.Ml. 4 d (the  maxi- 

mum rate  of  movement  by  the  total  pitch  d90ldt,  of  the  wobble  plate; 
the  wobble  plate  tilt  d9,ldt,  ; the  gradient  d 9,/dt  and  the  deflection 
d90Jdt).  of  the  total  pitch  of  the  tail  propeller).  Every  estab- 
lished rate  can  be  expressed  as  a certain  percentage  of  limiting 
rate . 

An  example  of  the  input  data  for  the  n-segment  is  given  in  Table 
1.  The  horizontal  line  in  the  table  of  data  means  that  the  given  para- 
meter is  a resultant  (passive)  one,  which  provides  only  the  informa- 
tion concerning  the  course  of  movement  and  which  is  not  included 
into  the  boundary  conditions  or  limitations.  Such  resultant  (pas- 
sive) parameters,  which  from  the  physical  point  of  view  (e.  g.  the 
perception  of  sense  and  the  pilot's  working  performance)  are  not 
able  to  serve  as  constraints  or  as  final  conditions  of  the  flight 
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Table  1.  The  input  datum  table.  For  the 
helicopter  exemplary  acceleration  at  the 
z0=3  m permanent  height  from  0-30  m/s  with 

the  0>-15°  impassible  gradient  as  well  as  a 
further  climb  on  the  z =15  m height  at  the 
80  km/h  speed. 
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Key:  (1)  Flight  Trajectory  Points;  (2) 

Location  - Coordinate,  Linear;  (3)  Velocities 
(Linear);  (4)  Location  Coordinate,  Angular; 
(5)  Boundary  Steering  Rate;  (6)  The. program 
of  Passage  from  the  n-segment;  (7)  The 
Program  of  Passage  from  the  n+1  segment. 

right  part  of  Table  1 on  following  page. 
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period,  are  the  values  ax. 


a , a which  belong  to  the  group  of  opera- 

y z 


tional  parameters;  w ,0)^, 


e , £ belong  to  this  group  as  well 

y z 


because  of  the  impossibility  of  physically  evaluating  their  magnitude 
during  the  flight.  however  the  above-mentioned  values  can  be  in- 
cluded in  the  group  of  data  in  the  case  of  piloting  by  means  of  au- 
tomatic pilot,  which  can  also  react  to  the  changes  in  derivative 
displacements,  or  in  the  case  of  numerical  analysis  of  a maneuver. 

At  times  a^  can  serve  as  a boundary  parameter  in  a "saddle  pressure" 
or  accelerometer  maneuver,  which  can  happen  in  the  case  of  sharp 
turn  or  pull-up. 


The  given  scheme  is  an  approximate  outline  of  the  system  of  the 
input  block  and  should  be  treated  as  a reference  tool  with  flexible 
application  to  the  tasks  dictated  by  the  intended  maneuver. 

2.1.2.  The  PN  Procedure.  Computation  of  Necessary  Power  During 
the  Unsteady  Helicopter  Plight 


The  acceleration  values  are  taken  from  the  FRn  program.  The 

helicopter  elements'  loading  is  taken  from  the  RS  procedure;  the 

KP  -P- 

motion  parameters  at  point  i are  taken  from  th^FN  procedure  is 
intended  to  calculate  the  essential  power  which  is  being  received 
by  the  helicopter.  Setting  up  the  quasi-stationary  system  which 
is  based  on  the  principle  that  the  induced  velocity  changes  are 
considerably  faster  than  the  changes  in  the  parameters  of  helicopter 
motion,  the  necessary  power  can  be  computed  in  the  following  way: 


the  resultant  thrust  of  the  helicopter  main  rotor: 


J 


r-Vrl+ij+Tl 

where  T , T , T are  components  of  the  thrust  in  the  helicopter 
x y z 

system. 

The  thrust  coefficient: 
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f 


Cf  = 0U*R5,nkbo 


where:  Q 

is 

the 

air  density; 

U 

is 

the 

speed  of  the  end  of  the 

rotor  blade; 

R 

is 

the 

radius 

k 

is 

the 

number  of  blades; 

bo 

is 

the 

chord  of  the  blade  root 

for  r=0 

The  induced  velocity 


where : 


v =.  \l  yi+  vl  + y\ 
v°  “ |/t«7TV 


is  the  flight  speed; 

is  the  speed  induced  in  hovering; 


x - coefficient  of  the  boundary  losses. 


For  low  flight  speeds,  for  — < 1.  in  order  to  allow  for  the 

»o 

non-uniformity  of  the  induced  velocity  field,  it  is  necessary  to 
assume  that 


*. 


and  for 


*1 


1 
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The  flow  velocity  through  the  main  rotor: 

~ Ksin  t+p  ^ the  dimensionless  flow  velocity  *= 
where  x is  an  angle  of  the  inflow  on  the  plane  of  the  main  rotor. 


Tangential  flow  around  the  main  rotor. 

V,p=VcosT  (the  dimensionless  velocity  of  the  flow  M,p=U). 

The  ground-proximity  influence  for  slight  tilting  <J>  of  the  main 
rotor  in  relation  to  the  Earth  can  be  evaluated  by  means  of 
factor  (7). 

K 1+  _ <?5.73  kb0ii  U2(y<p+v) 

~rr 


64TR3(hw  + z)2 


H* 


where:  bg  ^ is  the  chord  of  the  main  rotor  blade  for  r=0.7 

hw  is  the  height  of  the  main  rotor  hub  over  the  helicopter 
central  propeller; 

is  the  height  of  the  helicopter  central  propeller  over 

the  ground. 

The  ventilating  power^in  the  land 


The  profile  power  of  the  main  rotor  (as  a result  of  this  rotor's 
rotation) 


proximity/ 
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pwi 


1 0U3Rkbos(tA+  yM?  '») 


as  a result  of  the  tangential  stream 


Pwi  - -^-(0,425^t U +iir'ii  + 3/,’*)«l/JW:&o  K 

O 


where,  the  blade  profile  drag  coefficient  is: 


, rl+  w'crL  1 


Increase  in  power  due  to  the  compressibility  effect  AP 


Ecceeding  the  Ma  number  at  the  tip  of  the  blade 


AM  «=  — / -Mk 


If  AM  < 0 then  AP,  ■»  0 

If  AM  > 0 then  AP,  - 0,11 5/4  AM1! 5,73 kb. 


-'-xA- 


end  losses 


where  cr 


2 T 


tU^nA* 


The  blade  geometric  parameters 
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“ J -j—r*~xdr  — the  blade  contour  parameter 

<?.)  ° 

*1 

dr — the  blade  torsion  parameter 


Increase  In  power  as  a result  of  stream  separation  APq 

Having  the  lower  and  upper  Hy  separation  limits  one  can 
evaluate 

hp0  m JLrMt-  r 

0 Hv-Htwi 

Consequently,  the  total  power  consumed  by  the  main  rotor  is: 

The  power  consumed  by  the  tail  propeller 
The  tail  propeller  is:  T,  - -jr j- 

Ut  lt 

where:  l is  the  distance  between  the  tail  propeller  and  main 

rotor  shaft. 

The  velocity  of  the  inflow  on  the  tail  propeller  is  V^Vsini^, 
where  t is  the  inflow  angle  on  the  tall  propeller  plane  (the  blade 
tip).  The  flow  velocity  through  the  tail  propeller  is: 


„ v, - »*.„  - £ ♦ y'%  ♦ 


For  the  inflow  whose  direction  is  different  than  the  tail  propeller 
induced  velocity 


I Bfl 


* 


' 


2 )y,--v>+v,l  ( where:  is  deter- 

mined from  Figure  5) 

T 

where  p - — ~ and  consequently 

^/2tf.  where  Rfc  is  the  tail 
propeller  radius. 

The  tail  propeller  Induced  power 

Pu  =■  r,  v, 

The  tail  propeller  thrust  coefficient: 

87* 

c*'  “ tu}R,\nk,b91; 

where:  Ufc,  kfc , bQ  _ - respectively  are:  thd  blade  tip  speed; 

number  of  blades,  the  chord  of  the  tail  propeller  blade  at  0.7  R 

u 

The  profile  drag  coefficient  of  the  tail  propeller  blade  is: 

«>■«•  ..  I 

L *Jr+/*»(2fi,  + /|, /Jf)  J 

where: 

"i 

The  profile  power 

^ *•  *••*» *•  y hij 

The  power  consumed  by  the  tail  propeller. 


Figure  5.  Evaluation 
of  the  induced  velo- 
city of  the  tail  pro- 
peller in  the  area  of 
the  vortex  ring. 
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p,  * p„+p, 


The  total  power  necessary  for'  the  helicopter  flight: 

#>  = (/>„+/>,){ 

where  ^ is  a coefficient  of  the  power  utilization  with  consideration 
of  losses  connected  with  propulsion  of  assemblies. 

The  above-mentioned  formulas,  although  approximate,  will  be  used 
in  an  exact  model  for  the  evaluation  of  power  increases  arising  from 
small  changes  in  motion  parameters.  Therefore,  we  should  expect 
precision  an  order  of  magnitude  greater  than  with  the  evaluation 
of  0 states. 

2.1.3.  Procedure  N.  Computing  the  lift  of  the  main  rotor  T^ 
and  separation  ceilings  and 

Approximation  of  the  critical  lift  curve  CwU=f(y)by  analytical 
dependence : 

CwV  = 0,02  [29,8 - j/26j2* -20yr^2376)r+  (0,5 M0)  (0, 1 7 - 0,38//r)  + 

+ 0,0  57,3  - 7)  (0,0005  - 0,00235/<r)- 0,008/1]  (1  - /!)  - 0,005r//rw  57,3  + 

+ («*„- 14)  [(0,01 8 - 0, 024 A/0)  (1 -2//r)  + 0,004yr] 

C„l  “ CwV  — 0,1  +0,2yr 

if  pt>0.5  then: 

CwL  - C*,  - 5 (0,7-„r)  [0,02  (29,8-^26,2*-  13,6J)]-(0,-53/o)  0,02  + 

+Ort  57,3  - 7)  (0,0005  + 0,00235  0,5) -0,008/!  [(1  -O~0,005t  0,5  574+ 

+ («,»- 14)  0,004  0,5] 
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the  air  density  at  the  separation  ceiling 
the  air  density  at  the  separation  ceiling 

Ceilings  of  the  flow  separation: 

Hu  - 20  (0,131  -0,0004  fOjJ-<rl,)/(0,l 31  -0,0004  toll+e„) 

Hl  - 20  (0, 1 3 1 - 0,0004  t0il~  ql)I(0,  1 3 1 - 0,0004  toll  + ql) 

The  critical  load  carrying  ability  at  altitude  H 

1 

*'2  Q ^o,?  RU* 


Qv 


Ql 


2 Tw  

CmV  kb0.i  RU1 

CwV 
Qv  r' 

'-wL 


Symbols : 


M0  - 


U 

a* 

a*  — 

3.0- 

«*»,— 

foil 


the  Ma  speed  of  the  blade  tip  while  hovering; 

sonic  speed  altitude  H; 
geometrical  torsion 

flow  separation  angle  at  the  blade  tip; 
temperature  near  ground  for  given  type  of  climate. 


2.1.4.  The  TET  Procedure.  The  location  of  the  control  system 
of  the  main  rotor. 


As  a result  of  solving  equations  of  equilibrium,  the  location  of 
the  steering  system  of  the  main  rotor  is  calculated. 


Assumptions : 

- correction  of  the  reverse  flow  only  for  the  thrust;  any  blade  geo- 
metry is  included  by  using  the  kn  and  tn  coefficients;  flexibility 
of  blade  attachment  is  included,  making  it  possible  to  compute 
unarticulated  rotors  by  introducing  eauivalent  rigidity  in  a flapping 


J 


t 
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hinge  (according  to  Figure  6). 


- reverse  flow  for  blade  hinge  and  antitorque  is  not  included 

- simplified  geometry  for  the  induced  velocity  field  is  assumed 
(according  to  Glauert); 


- it  is  assumed  that  flow  around  a blade  element  is  perpendicular 
to  the  blade  axis; 

- constancy  of  the  gradient  of  lift  is  assumed  only  the  first 

dot 

harmonic  of  blades  motion  is  considered; 


For  a given  value  of  thrust  T,  rotor  inclination  x,  rolling  moment 
Mx,  pitching  moment  My,  flow  V^,  flow  V,p,  helicopter  angular  velocities 
<jjx  and  which  were  computed  as  a result  of  solving  equilibrium 

equations  RS  in  a given  flight  state,  the  following  items  are  computed: 

angle  of  blade  total  pitch 

# + *,+  \ r^*,+(/,  + 24)>l0 

®o  : 

*»+  2 fot i 


position  of  plane  of  constant  angles  of  incidence 

90—2ki  — iI  A0)+  -£■- 

- + 

t*~  f fo  ti 

V “ »•/ 

b foh*o+t4vK-Jjfz  / 4ci^  8w,  \ 1 

/.*;  rf<7  ' ri(ri+  j,!) 
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} 


where  y - 4/y*,  w = v/U,  K = 4/<r/[3  (l,2/!T+^r)]  y*  = MpIqR4  5,73  b0,  /p—  moment  of  inertia 
of  blade  in  relation  to  flapping  hinge. 


rotor  coning  angle: 


where  C»*  - HkH/eU2R25,73  b0. 


and  a^  are  computed  by  the  method  of  successive  approxima- 
( tions;  assuming  the  first  approximation  for  a^. 

a\0)  = Hj  — (2C*  + t2  It) 

U 


we  compute  ft*  So"  for  the  assumed  aj^  and  then  and  a^1^  until 
achieving  an  accuracy  of 

\a[J)-a<J-t}\  *5  0,00873 


\T' 


Figure  6.  Scheme  for  replacing  non- 
articulated  rotor  and  flexible  blades 
with  articulated  rotor  and  rigid  blades 
and  with  equivalent  flexibility  k in 

flapping  hinge:  1-actual  blade  of  non- 
articulated  rotor;  2 - substitute 
blade.  (a)  substitute  blade;  (b) 
rotation  axis;  (c)  actual  blade  of 
non-articulated  rotor. 
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f 

# 

which  usually  happens  in  the  2nd  or  3rd  approximation. 

the  flow  components  through  the  rotor  on  the  plane  of  constant  angles 
of  incidence: 

V.  - VA  cos  a,  + yT  sin  <»i 

V,  = — VA  sin  a,  + VT  cos  a, 


dimensionless  values 


A0 


R o = 


V. 

a)R 

Vj_ 

a)R 


wobble  plate  location  in  relation  to  the  shaft 

9„  = [(a,,  - €>  + y,)  Z>,  - (4>  - y,  - *„)  Z>2]/(/>?  + 

where : 

D1-kDl 


Dl  + kD1  . 

- “V-TF  ' "■  -i+P 


D^,  - the  control  system  coefficients  (data  from  control  system 

design) ; 

k - the  kinematic  connection  of  oscillations  and  turns. 


Location  of  the  wobble  plate  in  relation  to  the  plane  of  the  main  rotor 
blade  tips. 
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aly  = 


a | +ibt 
l+fc1 2 


bi~kat 

\+i 


inclination  and  bending  of  the  thrust  vector  of  the  main  rotor. 

T. 

V,  = arctg  jr- 
T 

V,  = arctg-^ 


2.1.5.  The  RS  Procedure.  The  Helicopter 
Equilibrium  Computation 

The  most  important  problem  as  well  as  initial  action  is  to  bring 
the  motion  equations  into  agreement  with  the  physical  equivalent  of 
the  phenomenon  which  is  being  modeled.  The  equation  system  has  to 
consider  both  the  evaluated  motion  parameters  and  the  coefficients 
which  connect  them  depending  on  the  helicopter  design  parameters  as 
well  as  on  its  system  characteristics.  The  motion  equations'  system 
in  space  has  to  consider  the  sum  of  forces  and  moments  in  relation 
to  three  axes  and  it  should  also  consider  the  equation  of  motion  (or 
energy)  of  Internal  degrees  of  freedom,  primarily  of  the  system  of  the 
main  rotor.  In  most  cases  In  modern  helicopters  the  additional  motion 
equation  has  a lower  value,  since  (contrary  to  airplanes)  in  case  of 
employment  of  amplifiers  in  the  rotor-to-controls  direction  the  system 
is  self-braking,  which  makes  superfluous  the  evaluation  of  the  dynamics 
of  motion  with  a "release  rudder".  To  evaluate  the  above-mentioned 
dynamics  an  addition  equation  of  the  control  system  was  necessary^ 
in  the  general  case  the  motion  equations  must  include: 

1)  the  motion  equation  of  a helicopter  as  a rigid  body  in  a system 

with  an  axis  connected  with  the  helicopter  (the  Eulerian  equation), 
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2)  equations  of  motions  of  control  systems; 

3)  the  motion  equation  (energy,  power  balance)  of  the  inertial  system 
of  the  main  rotor; 

4)  possible  equations  of  deformation  of  the  helicopter 
The  equation  of  the  motion  of  the  helicopter  element  <5 

m 

SF  - Sm  dv 
dt 


Completing  the  integration  we  will  obtain  the  force  equation: 


and  from  here  we  will  obtain  3 scalar  force  equations 


&h  = r vSm 

The  elementary  moment  of  motion 


and  from  here  3 scalar  moment  equations. 


For  the  co-ordinate  system  according  to  Figure  7 one  can  obtain  6 
motion  equations. 


Figure  7.  The  accepted  co-ordinate 
system.  Key:  (1)  Center  of 
gravity. 
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Beginning  with  a system  of  axes  connected  with  the  helicopter  at 
the  center  of  gravity,  we  will  obtain  the  following  equation  system: 


m (a,  + V,  to,  - F,  <o,)  + mg  sin  0 = X 
m (a,+  F,  to,  -to,  V,)  — mg  cos  0 sin  <P  = Y 
m (a,  + F,  to,  - F,  to ,)  — mg  cos  9 cos  4>  = Z 
I,  + (/,  - lf)  to,  to,  - /,/«,  - to,  to.)  + /„(wj  - to*)  - /,,(e,  + to,  to,)  = L 
1.)  to,  to,  - /„(e,  - to,  to,)  + /„(to?  - col)  - /„(e,  + to,  to,)  = M 

, 

/,  e,  + (/,  - /,)  to,  to,  - /„(£,  - to,  to,)  + /„K2  - to,2)  - /„(£, + a>,  to,)  « AT 

The  aero-dynamic  forces  and  moments  located  on  the  right  side 
of  the  equations  are  non-linear  functions  of  the  flight  speed,  wobble 
plate  deflection,  main  rotor  total  pitch,  tail  rotor  pitch  and  gust 
speed. 

The  angular  velocities  in  relation  to  the  axes  connected  with 

the  helicopter  and  designated  in  equations  as  oo  , to  , u>  are  also 

x y z 

non-linear  functions  of  the  Euler  angles  *,©,  v 

1 • 

w,  - 4- y sine 

ft),  - 0 cos  <►  - ¥ COS  0 tin  * 
ft),  ■«  IP  cos  9 cos  sin  # 

The  rotor  equation  of  moments  (power) 

P'  + Pt  + P.'+P,  = max  Vx+GV,  + l(oa>+P, 

where:  Pp  and  pj^  are  the  profile  and  Induced  power  of  the  main  rotor j 

Psz  ls  P°wer  lost  to  parasitic  drag; 

-- 

PQ  is  power  consumed  for  climb  or  descent; 


27 


V is  vertical  speed  (negative  for  descent); 

I is  inertial  moment  of  the  inertial  system  of  the  main  rotor 
reduced  to  the  rotor  axis  and  revolutions; 

o>  is  the  rotor  angular  velocity; 

is  flight  speed  (horizontal); 

ax  is  the  helicopter's  acceleration; 

is  power  supplied  from  propulsion. 

The  general  solution  of  the  equation  system  can  be  achieved  by  the 
numerical  method.  It  is  applicable  in  cases  when  the  decisive  elements 
of  motion  remain  unknown  and  simplifying  assumptions  cannot  be  applied. 
Quite  often  it  is  necessary  to  analyze  significant  motion  disturbances 
(e.g.  when  evaluating  maneuvers),  where  non-linear  dependences  and 
coupling  of  equations  cannot  be  forseen.  The  numerical  solution  process 
through  the  operations  performed  on  numbers  instead  of  formula  conver- 
sion) of  the  helicopter  transient  equilibrium  is  an  approximate  process. 

The  "step  by  step"  method  is  applied  (see  the  KR  procedure). 

Starting  from  the  V , V , V , u>  , u>  , u>  initial  values  of  the  velo- 

x y z x y z 

city  as  well  as  the  rudders'  locations  and  the  rotor  thrust  components 

it  is  necessary  to  compute  the  forces,  X,  Y,  Z,  M , M , and  after  that 

y z 

the  values  of  acceleration  ax,  a^,  ag,  ex,  e , ez.  Then  assuming  that 
the  above-mentioned  accelerations  during  time  period  At  are  constant,- 
the  new  values  of  velocity  obtained  by  the  end  of  time  period  At  are 
computed.  Accepting  the  new  speed  values  as  the  output  ones  for 
computations  during  the  next  time  period  At  one  can  compute  in  the 

step-by-step  way  the  entire  pattern  of  motion  while  considering  the 

♦ - 
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Figure  8.  The  helicopter  orientation  in  space.  Key: 
(1)  flight  trajectory;  (2)  Perpendicular;  (3)  System 
connected  with  the  ground. 


Figure  9.  The  motion  parameters  computed  in  the  KR  phase. 


influence  of  change  of  forces  and  without  limiting  oneself  to  the 
case  of  the  linear  dependences  of  forces  on  the  steady  motion 
disturbances  (which  is  very  important  with  significant  flight  devi- 
ations or  in  flight  conditions  with  flow  separation  or  within 
the  compressibility  influence).  Changing  time  period  At,  one  can 
achieve  more  precise  computation  results. 

2.1.6.  The  KR  Procedure.  (The  Solution  Step  for  the 
Time  Interval  At  Between  t-^  and  t^+^ 

The'  segment  of  the  flight  trajectory  with  a uniform  transition 
technique  points  n-n+1  should  be  divided  into  a series  of  time 
intervals  for  precisely  determining  the  pattern  of  variation  in 
parameters  of  motion  between  points  n and  n+1. 

This  procedure  deals  with  the  problem  of  these  parameters’ 
behavior  within  the  n-n+1  segment  in  order  to  control  the  impulse 
of  the  thrust  components  ATx,  ATy,  ATz,  AT^  of  the  rotor  which 
are  brought  into  being  through  the  displacement  of  the  control 
system  by  AS,,  AS,,  A90,  A3,.  The  accelerations  of  the  helicopter,  which 
is  assumed  to  be  constant  in  the  time  At,  correspond  to  increments 
in  thrust.  These  accelerations  are  equal  to  the  accelerations  at 
point  i or,  in  the  case  of  the  iterative  repetition  of  computations 
within  an  interval  of  one  step,  to  the  mean  value  computed  with 
the  following  formula: 


<axi)+axiil)/2 

where  j is  a succeeding  step  iteration 

At  the  beginning  and  end  of  the  i-step  the  following  parameters 
are  analyzed:  the  location,  speed,  and  acceleration  of  the  helicop- 
ter; the  values  of  the  aerodynamic  and  mass  load  of  the  helicopter 
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components  as  well  as  the  power  consumption  through  the  rotor  system. 
The  angular  location  (♦,©,¥')  is  necessary  in  order  to  analyze  the 
wobble  plate  location  in  relation  to  the  rotor  shaft  for  the  purpose 
of  the  determining  steering  reserve  in  the  following  step  and  the 
possibility  of  bringing  about  a proper  impulse.  The  wobble  plate 
possesses  a small  range  of  angular  deflections  (in  relation  to  the 
shaft),  and  consequently  the  helicopter  location  in  space  influences 
significantly  the  control  over  the  value  and  inclination  of  the  thrust 
vector. 

For  illustration  certain  parameters  of  the  helicopter  equilibrium 
computed  in  the  KR  phase  are  shown  (See  Figure  9). 

For  small  ^t  and  the  average  values  of  the  accelerations  ax,  a , 
x^,  ex,  £y,  ez,  the  parameters  of  the  end  of  step  i+1  can  be  expressed 
through  the  parameters  of  the  beginning  of  the  step: 


xi+t  — xi+ 

At2 

y 1+1  = yi+  Y,i  At  + af  2”- 
A t2 

zi+i  = * i + Ki  At+a2  2 
At2 

*i+i  = <t>i  + (Dx,At  + ex  2 
At2 

0(+,  = 0,  + vt,,  At  + c,  2 
At2 

Ki+ 1 = V„i  + axAt 
= y,i  + arAt 
^«i+i  “ y„  + a.  At 
«**/+!  " (oMl+eMAt 
*>H* 1 - (o„+e,At 
».«♦  i - a»„  + «.  At 


displacements  at 
point  i+1 


velocities  at  point 
i+1 
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2.1.7.  The  D Procedure  (derivatives,  control  gradients) 


0.  T,  = /(->„,) 

As  a result  of  the  solution  of  the  helicopter  equilibrium  equations 
by  means  of  the  RS  procedure  as  well  as  by  means  of  the  changes  in  the 
motion  parameters  evaluated  in  the  KR  procedure,  the  location  and  load- 
ing of  the  main  rotor  and  the  flow  around  it  are  determined. 

V V Tz  " load; 

VA,  VT,  VB  - flow; 

oj  z - the  angular  velocities  of  helicopter  roll, 
pitch,  and  yaw  in  relation  to  the  center. 

There  are  no  methods  yet  for  the  precise  computation  of  the  depend- 
ences T . T , T “/(*,.  9*.  #o)  for  the  given  parameters  of  flow  around 
x y z 

and  motion  of  helicopter  in  space,  which  could  serve  for  the  determi- 
nation of  the  derivatives  <fr,1  Jr'  , etc.  in  order  to  evaluate  the  thrust 

dt>u  Wtt 

impulse  caused  by  the  wobble  place  motion  with  values 

However,  for  the  given  values  of  load,  flow  and  angular  velocity 

of  the  rotor  in  space  it  is  possible  to  find  relatively  precisely  the 

increments  A80.  Aft,.  A9,  for  the  given  AT,,  AT,,  AT,,  or  to  solve  an  inverse 

problem.  It  allows  the  following  procedure:  the  necessary  power  and 

the  positions  of  the  control  systems  at  a given  level  of  the  rotor 

load  are  computed.  Then,  for  a given  increment,  e.g.  AT  = 100  kG  (or 

2 z 

acceleration  impulse,  e.g.  Aa  1 m/s  ) new  powers  and  positions  of  the 

z 

control  systems  are  computed  and  from  this  we  obtain  such  required 
derivatives  as: 


4,  - \P ■ 

Afl, 


where 


” P«T.+*T.)  — P, 


■«r> 


t 


\ 

I 

f 

AT*  s 

dr  - where,  for  instance  A7,  = 100  kG  is  assumed. 

These  derivatives  will  be  used  both  in  the  evaluation  of  the  margins 
of  control  and  in  the  selection  of  the  permissible  control  impulse 
in  the  i-step. 


2.1.8.  The  PR  Program.  Flight  on  n-f light 
Traj  ectory 

This  program  represents  a mathematical  model  for  the  method  of 
accomplishing  of  a segmental  task.  It  is  constructed  to  correspond 
physically  to  an  actual  pilot-helicopter  system  and  piloting  techniques  \ 

used.  In  the  program  of  flight  for  segment  n there  must  be  a coded 
principle  of  this  flight  as  well  as  the  goal  realization  method,  the 
so-called  final  parameters  of  the  segment.  Formulating  then,  according  ‘ 

to  the  physical  sense  and  realization  technique  of  a proper  maneuver, 
the  mathematical  model  of  its  realization,  attention  should  be  drawn 
to  the  following  problems: 


1)  energy  distribution  (allocation  of  the  available  power  to  climbing, 

| 

horizontal  acceleration  and,  if  necessary,  to  acceleration  of  the  rotor, 
etc.)  which  in  reality  is  accomplished  intuitively  by  the  pilot,  who 
sets  the  values  of  the  parts  arbitrarily; 

| 

2)  connecting  the  maneuver  segments  (so-called  transitional  states) 
between  two  different  sequential  maneuvers; 

3)  forseeing  the  helicopter  motion,  that  is,  the  computation  and  extra- 
polation parameters  of  the  current  moment  in  order  to  evaluate  the 
probability  of  the  reaching  goal  and  to  correct  the  actual  control 
Impulses  in  such  a way  that  the  given  parameters  of  the  end  of  the 
segment  could  be  obtained  with  a permissible  tolerance; 

4)  the  tolerance  for  achieving  the  final  parameters  of  motion; 


5)  selection  of  the  helicopter  control  impulse  during  step  i. 

The  co-ordinate  cuboid  (Figure  10),  at  whose  diagonal  ends  the 
boundary  conditions  of  helicopter  motion  are  determined,  is  a unified 
element  of  the  flight  profile  which  makes  it  possible  to  realize 
virtually  any  task,  which  can  be  divided  into  segments  characterized 
by  similar  pilotage  techniques.  There  is  always  a certain  option  in 
switching  from  point  n to  point  n+1  but  in  the  algorithmized  procedure 
there  must  be  an  unequivocally  coded  method  of  reaching  the  required 
motion  parameters,  taking  as  a basis  the  parameters  given  at  the  n- 
point.  This  option  is  limited  by  the  boundary  conditions,  the  assumed 
rate  of  control,  and  the  accepted  principle  of  the  flight  for  segment 
n.  If  it  is  difficult  to  achieve  the  interconnected  parameters  of 
terminal  point  n+1  in  one  maneuver  computation,  the  computation  should 
be  repeated  and  the  solution  obtained  by  the  iterative  process.  Cases 
of  this  kind  occur  very  often  when  the  flight  task  is  being  optimized. 

2.1.9.  The  F System  (prediction). 

The  Model  of  the  Pilot’s 
Action 

The  system  is  to  evaluate  the 
motion  parameters  in  the  immediate 
future  and  to' provide  information  on 
whether  and  in  what  way  to  change  the 
helicopter  control  functions,  based 
on  the  anticipated  values  of  the 
motion  parameters,  in  order  to  obtain 
the  given  characteristic  values  of  the  function  of  the  goal  (the  para- 
meters at  point  n+l).  It  is  a physical  equivalent  of  the  pilot’s 
earliest  reaction  according  to  the  anticipated  motion  in  order  to 
obtain  a definite  flight  condition.  This  system,  by  extrapolating 
the  current  parameters  of  motion,  evaluates  the  effect  of  the  control 


r 


Figure  10.  Option  flight  paths 
from  point  n to  point  n+1. 


system  action.  Basing  itself  on  the  evaluation  of  the  situation  in 
the  step  which  is  under  linear  consideration  or  onA approximative 
function,  this  system  predicts  the  further  pattern  of  motion.  Thus, 
by  correcting  the  c on trol^e very  step  and  advancing  correspondingly, 
it  is  possible  to  reach  the  specific  motion  parameters  at  the  moment 
when  the  maneuver  n is  over. 

In  the  controlled  objects,  depending  on  efficiency  of  control 
and  the  actual  parameters  of  motion  one  can  separate  the  limiting 
time  interval  where  there  is  a possible  required  change  in  parameters 
of  the  motion  when  the  limiting  possibilities  of  control  are  used. 

This  value  is  of  a significant  importance  in  goal-oriented  control  when 
the  pilot's  reaction  is  based  on  an  anticipated  flight  state  and  not 
on  a current  set  of  motion  parameters. 


2.2.  Theoretical  Bases  for  Determining 
limiting  time  At 

In  the  general  case  it  is  the  time  necessary  for  a change  in  the 
vector  of  velocity  with  value  A^  (Figure  11)  where  is  the  difference 
between  the  velocity  of  the  goal  target  at  point  n-^n  and  the  velocity 
at  point  i-V.  with  the  assumed  resolution  of  the  maximum  acceleration 

1 'A 

in  segment  i-n.  (a)o 

AP  0 

\ 

The  resultant  motion  of  the  helicopter  in  a maneuver  is  a super- 
positioning of  displacements  in  relation  to  its  axis.  For  every  degree 
of  freedom  one  can  anticipate  the  helicopter  shifts  in  the  computational 
program  and  calculate  the  limiting  time. 

For  example,  the  physical  picture  of  changes  in  the  boundary  zone 
which  determines  the  prediction  time  necessary  for  further  prediction 
of  motion  is  given  in  the  example  of  the  change  in  the  vector  of 
velocity  at  point  n at  90°  during  accelerated  flight  from  point  n-1 

to  point  n,  starting  from  velocity  V = 0 (Figure  13). 
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Figure  11.  Change  in  the 
vector  of^velocity  AV  between 
velocity  Vn  and  the  velocity 

at  a current  instant. 


Figure  12.  Elementary  maneuvers  as  a basis  for 
determining  limiting  time  At  . 


Figure  13.  Change  in  the  region  and  limiting  path  of  the 
change  in  velocity  with  a change  in  flight  velocity  along 
the  path:  1 - limiting  path;  2 - penetration  boundary 
zone;  (3)  boundary  zone  radius. 
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Figure  14.  Figure  15. 

Figure  14.  Evaluation  of  the  range  of  permissible  change  in  the 
flight  path.  Key:  (1)  inaccessible;region;  (2)  limiting  path; 
(3)  of  goal. 

Figure  15.  Determining  At  as  a transposition  of  the  solution 
of  the  flight  path  from  Figure  14. 


The  necessary  limiting  time  can  be  computed  v.ith  the  following  il- 
lustration of  a change  in  the  direction  of  flight  (Figures  14,  15). 

There  is  a similar  representation  of  the  value  of  acceleration  or 
braking  of  linear  or  angular  movements  (Figure  16).  The  scheme  of 
action  of  system  F can  be  illustrated  in  the  example  of  a helicopter 
vertical  climb  to  the  preset  altitude  Zn=H  (Figure  17).  In  this  case 
x1=y^=z^=X2=y 2=0 . Continuous  computation  of  the  possibility  of  braking 
the  motion  makes  it  possible  at  the  proper  time  (point  A)  to  decide 
about  braking  the  motion  and  at  the  point  B to  decide  about  the  reaching 
the  given  altitude  of  climb,  which  in  this  simple  maneuver  was  the  goal 
and  the  boundary  condition  of  the  end  of  segment  n+1. 

When  setting  the  F system  procedure  one  should  consider  the  fol- 
lowing elements:  1 - evaluation  of  the  difference  of  the  current  and 
goal  parameters  and  the  instantaneous  limiting  margins  of  control; 

2 - the  principle  of  achievement  of  the  goal  function  (the  condition 
for  correction  of  the  control  steering  system);  3 - the  limiting 


; 

I 


I 


t 


37 


' 


. 

1 


Figure  16.  Evaluation  of  At  using  as  an  example  the 

horizontal  flight  acceleration. 

Key:  Maximum  possible  acceleration. 


possibility  of  a change  in  motion  in  the  immediate  future(for  estimating 
of  the  necessary  time  for  prediction  of  Atgr). 

2.3.  Distribution  of  Excess  Available 
Engine  Power 

Program  PRn  often  involves  the  problem  of  how  to  allot  engine  power 
in  order  to  obtain  the  conditions  given  at  the  end  of  a segment  for 
the  maneuver  under  consideration.  In  the  general  case  the  distribution 
of  power  in  three-dimensional  motion  can  be  in  5 directions:  1)  for 
climbing;  2)  for  horizontal  acceleration;  3)  for  lateral  acceleration; 

4)  for  accelerating  the  rotor;  5)  for  the  tail  propeller  in  case  of 
rotation  around  the  vertical  axis.  The  necessary  power  changes  cause 
changes  in  the  flow  around  the  rotor  and  in  the  rotor  loads  during  an 
instantaneous  state  of  helicopter  equilibrium  during  a maneuver.  The 
kinetic  energy  of  the  rotor  Inertial  system  can  be  considered  as  a 
part  of  the  available  power  because  of  the  specific  nature  of  its 
source . 

The  shares  of  the  excess  power  au  in  each  of  the  5 above-mentioned 
directions  can  be  presented  as  projections  of  a unit  vector  of  control 
in  a five-dimensional  system. 
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Figure  17.  Function  of  system  F in  single-parameter 
helicopter  motion. 

Key:  (1)  Decision  to  brake  in  motion  taken  according 

to  system  F;  (2)  Goal  - climbing  ceiling,  Limiting  path. 

For  the  first  approximation  one  can  often  apply  an  approximative 
analytical  formula  (e.g.  in  [7],  [10])  where,  assuming  average  values 
for  the  flight  phase  parameters  and  solving  the  system  of  equations 
for  uniformly  accelerated  motion,  one  can  compute  approximative  au 
values.  They  assure  achievement  of  the  parameters  required  in  a 
definite  time,  allowing  for  their  continuous  correction  in  each  1-step. 

2.4.  The  Problem  of  Transient  States 

In  order  to  carry  out  the  flight  task  it  is  necessary  in  the  PRn 
program  to  connect  individual  maneuvers  with  transient  phases.  The 
condition  is  that  the  final  parameters  of  the  preceding  segment  become 


Figure  18.  Figure  19. 

Figure  18.  Assumed  flight  path. 
Key:  (1)  assumed  flight  path. 

Figure  19.  Isolated  transitional 
phase. 

Key:  (1)  Transitional  phase;  (2) 

realized  path. 
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Figure  20.  The  transitional  phase 
Included  in  the  beginning  of  the 
next  segment  of  the  maneuver. 

Key:  (1)  Transitional  phase. 

the  initial  conditions  of  the  following  segment.  In  this  case  there  are 
2 ways  possible:  1)  employment  of  an  isolated  transient  phase  during 
which  the  above-mentioned  condition  of  the  conformabllity  of  the  final 
parameters  of  the  preceeding  segment  with  the  initial  parameters  of 
the  following  segment  according  to  Figure  19  is  realized,  2)  inclusion 
of  a transient  phase  in  the  beginning  of  a subsequent  segment  (Figure 
20). 


The  second  way  is  more  often  in  use  because  of  the  greater  ef- 
ficiency of  maneuvers  (the  maneuver  becomes  more  compact).  In  the 
phase  of  connection  of  procedures,  the  beginning  of  the  following 
segment  should  be  assumed  as  a function  of  the  goal,  and  the  factor 
adjusting  che  helicopter  flight  parameters  to  the  conditions  of  the 
following  segment  is  the  movements  of  the  flight  controls  at  a given 
control  rate. 

2.5.  The  Problem  of  Tolerance  in 
Achieving  the  Final  Parameters  of 
a Segment 

The  final  values  of  motion  must  be  achieved  with  a certain  toler- 
ance. The  tolerance  value  and  direction  should  also  be  connected  with 
the  kind  of  maneuver  assumed.  The  tolerance  value  (generally  consider- 
able)  isA individually  on  the  basis  of  experience,-  since  excessively 
tight  intervals  of  the  tolerance  value  prolong  the  computation  time 
(additional  iterations)  as  well  as  the  time  for  completing  the  maneuver. 

2.6.  The  Problem  of  Control  Impulse 
Selection 

In  the  control  program  (when  the  rotor  thrust  impulse  is  being 
chosen)  the  impulses  A*.  of  the  motion  of  the  helicopter  flight  controls 
are  evaluated  (Ax^  is  symbolically  understood  as  At*  A9X,  A3,). 

The  following  basic  constraints  must  not  be  exceeded  during  step  i: 

1)  available  power  Pr-Pn>0, 

2)  rotor  lift  (take-off)  Tkr~T>0, 

3)  design  constraints  of  each  channel  of  control  (applies 

to  collective  pitch,  cyclic  pitch  of  inclination  and  tilt  tail 
propeller  pitch); 

4)  the  requirements  of  the  flight  comfort,  e.g.  #>  -Iff 
(slope  of  passenger  helicopter  floor)  or  [4)<43g 
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(the  vertical  acceleration); 

5 ) dSJdt  < (d»,ldt)Ur.  of  the  permissible  rate  of  control. 


There  can  be  other  constraints  dictated  by  the  specifics  of  the 
maneuver  performed.  Prom  procedure  D one  can  obtain,  for  a given 
rotor  load  level  and  flight  state  computed  from  the  equilibrium 
conditions  (the  RS  procedure),  steering  derivatives  (linearized  for 
a low  impulse)  of  the  type  dT  * » etc. 

Hence,  the  maximum  control  margins  can  be  evaluated: 


" *,  ' 


K-T 

dj 


One  choses  the  minimum  margin  which  insures  that  the  given 

constraints  in  step  i are  not  exceeded.  On  obtaining  , one  can 

evaluate  the  value  of  the  increase  in  the  thrust  component: 


A 7)  m Jf 


i 


This  is  a first  approximation,  since  the  positions  of  the  wobble 

plate  are  precisely  determined  at  the  beginning  of  step  i and  the 

end  of  i+1.  Every  step  thus  begins  with  the  values  T , and 

x y , z 

of  determined  through  solution  the  inverse  problem.  It  provides 
a complete  picture  of  the  pattern  of  forces  and  motions  with  the 
flight  controls  as  a function  of  time  during  the  solution  of  the 
flight  problem.  In  practice  it  is  sufficient  to  consider  the  coupling 
o f 7j[ (Pot  Qf)  and  UK  »,),  while  the  coupling  with  and  the  depend- 
ence on  UK  K can  be  omitted. 


From  the  physical  point  of  view  this  means  that  the  rotor  inclina- 
tion (in  order^^o^^oi^ionent  T ) requires  the  correction  (increase) 
of  the  collective  pitch,  or,  when  the  collective  pitch  (in  order  to 
increase  force  T ) is  increased,  a correction  of  the  rotor  inclina- 
tion  (decrease  in  the  plate  inclination)  is  required  during  oblique 
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flight  in  order  to  maintain  equilibrium  conditions  along  the  x-axis 
(Figure  21). 


For  given  values  of  AT  and  AT  the 
evaluation  of  increments  A<?  n and 

o y 

can  be  obtained  through  the  solution  of 
system  of  2 equations  with  2 unknowns: 

Ar.-«4S"+-£48’+r-<4“) 
4r--^4a°+^49’+n(4’) 

Figure  21.  Diagram  of  load 
impulses  by  the  rotor, 

depending  on  the  action  of  where  the  increase  in  the  thrust  vector 
the  control  system  plane  . incllnatlon  ls 
x,  z. 

Act  «/(A9„AS0) 

- The  control  system  action  brings  a- 

bout  corresponding  increases  in  the 
thrust  components  in  relation  to  the  system  connected  with  the  heli- 
copter. Therefore,  in  order  to  analyze  the  flight  path  in  the  system 
connected  with  the  earth,  it  is  necessary  to  consider  the  change  in 
, the  reference  system. 

I 

3.  Exemplary  Illustration 

* 3.1.  A study  of  the  Operation  of  Elements 

of  the  Proposed  Method  to  Analyze  Non- 
Stationary  Helicopter  Motion  Taking  as  an 
Example  a Vertical  Maneuver  [9] 

• The  assigned  physical  behavior  of  the  phenomenon  is  as  follows: 

the  helicopter  climbs  vertically  from  point  Z^  to  point  Z2, 
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descends  from  point  Z 2 to  point  Z^,  and  climbs  to  point  Z^.  The 
points  designated  as  Zn(n=l,  2,  3,  4)  are  characteristic  points  of 
the  flight  path  (Figure  22). 


parameter  control  motion 
of  helicopter. 

Key:  Direction  of  mo- 

tion for  climbing. 


The  vertical  path  is  divided 

into  segments  with  equal  time  intervals 

At.  The  motion  is  solved  on  a "step- 

by-step"  basis.  The  beginning  of  a 

step  is  designated  z^  and  its  end  z^+-^. 

Due  to  the  discontinuity  of  the 

descreet  system,  at  the  beginning  and 

end  of  a step  one  can  distinguish 

the  left  side  of  section  z^  and  the 

right  side  of  section  z which  are 

P 

directed  along  the  direction  of  motion. 
The  solution  is  based  on  a group  of 
typical  procedures  combined  in  a common 
system.  In  one  computation  step  the 

Jp(i) 


right  side  of  step  i-z  , and  the 


transition  to  z£(  j_+i ) are  computed. 
The  principal  block  connecting  the 
entire  system  is  the  PRn  program,  in 


which  there  is  a coded  way  of  realizing 
helicopter  motion 

in  the  flight  segment  under  consideration.  Here  one  choses  a suitable 
increase  in  the  collective  pitch  of  the  wobble  plate  designated  as  A«?0. 
The  program  employs  auxiliary  procedures  and  systems,  such  as:  procedure 
Pn,l^  (computation  of  the  necessary  power  and  of  the  wobble  plate  posi- 
tion); procedure  KR  (computation  of  the  parameter  changes  during  one 
step);  system  F (the  area  penetration  in  the  immediate  future  and  pos- 
sibly the  correction  of  the  assumed  principle  of  control);  procedure  D 
(control  gradients),  etc.  The  solution  of  the  helicopter  motion  is 
shown  on  Figure  23. 
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Figure  23.  Pattern  of  changes  in  parameters  of 
motion,  control  function  and  loading  during  a 
simulated  single-parameter  helicopter  maneuver 
in  the  perpendicular. 

Key:  (1)  Rotor  thrust;  (2)  necessary  power; 

(3)  total  pitch;  (4)  acceleration;  (5)  speed;  (6) 
displacement;  (?)  time  partition;  (8)  engine  power; 
(9)  steering  rate. 


3.2.  Simulation  of  the  Helicopter 
Interrupted  Take-Off. 

The  method  usefulness  was  verified  by  modelling  a complex  heli- 
copter problem  - a take-off  with  engine  failure  [10]. 

The  simulation  model  of  controlled  flight,  as  applied  to  the 
analysis  of  the  helicopter  take-off,  includes  certain  simplifications 
namely : 


- only  two-dimensional  motion  in  the  plane  symmetry  x,  z is  ■ considered 


- the  helicopter  is  considered  as  a material  point  (the  equilibrium 
of  moments  in  relation  to  the  transverse  y axis  is  ignored); 

- it  is  assumed  that  the  helicopter  pitch  is  approximately  equal  to  the 
tnrust  vector; 

- in  view  of  the  low  analyzed  flight  speeds,  the  computational  algo- 
rithms omit  specifics  which  are  of  importance  for  high  flight  speeds; 

however,  the  following  must  be  considered: 

- the  ground  effect  as  a function  of  the  flight  speed  and  the  influ- 
ence of  the  inclination  of  the  rotor  plane  in  relation  to  the  ground; 

- optional  structural,  pilotage,  aerodynamic,  flight-comfort  and  other 
constraints  which  are  essential  from  the  point  of  view  of  helicopter 
take-off  technique; 

- the  tolerance  for  execution  of  maneuvers  in  the  take-off  phase  which 
correspond  to  the  pilot's  level  of  perception  in  estimating  actual 
flight  parameters. 

In  the  computational  program,  a special  control  program  has  been 
worked  out  for  every  typical  flight  segment;  it  evaluates  control 
impulses,  allows  for  constraints,  obeys  principles  of  flight  on  the 
given  segment  and  monitors  achievement  of  the  goal.  This  program  uses 
standard  procedures,  typical  for  all  flight  segments,  such  as:  the 
evaluation  of  the  helicopter  momentary  equilibrium  in  unsteady  flight; 
the  control  system  position  in  relation  to  the  helicopter. 

The  result  of  the  completed  simulation  process  is  the  profile 
of  the  helicopter  motion  parameters  which  is  shown  on  Figure  24 
where : 
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Figure  24.  Interrupted  take-off-landing  G = 4450  kG, 
H=0,  AW;  profile  of  changes  in  parameters  of  helicopter 
motion. 

Key:  (1)  Engine  failure;  (2)  time  for  pilot  reaction; 

(3)  pulling  up;  (4)  pulling  off;  (5)  flight  path. 
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V is  the  helicopter  horizontal  speed; 
w is  the  vertical  speed; 

a and  a are  the  horizontal  and  vertical  accelerations  of  the 

A Zi 

helicopter; 

is  the  collective  pitch  angle; 

<p  is  the  helicopter  pitch; 

U is  the  blade  tip  speed. 

Suggestions 

The  above  mentioned  employment  of  the  similar  analyses  can  be 
also  applied  to  other  helicopter  flight  problems.  The  information 
obtained  in  this  way  gives  a comprehensive  picture  of  the  pattern  of 
effects  during  the  task  being  performed,  with  easy  computation  of  i 

different  variants  of  the  flight. 


Errors  in  the  motion  evaluation  arising  from  the  adoption  of  an 
approximate  helicopter  model  are  of  no  special  significance,  since  in 
reality  the  repetition  of  maneuvers,  because  of  different  flight 
conditions,  pilotage  techniques,  etc.,  is  not  a strict  rule.  The 
spread  is  generally  large  - larger  than  the  error  of  the  method. 

In  view  of  the  ever-increasing  popularity  and  utilization  of  the 
method  of  representing  reality  through  the  fullest  possible  simulation, 
in  order  to  solve  the  problems  of  the  flight  simulation  (helicopter, 
glider,  airplane)  it  is  necessary  to  undertake  and  direct  studies  to 
determine  the  perception  of  the  flight  parameters  in  the  "man-machine" 
system. 


The  data,  mostly  in  the  field  of  pilot  psychophysiology  (e.g. 
reaction  time,  visual  estimation  of  flight  altitude,  etc)  must  be 


specified  numerically,  since  they  represent  the  input  data  in  the 
computational  program. 
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Summary 


A method  for  digital  simulation  vs •« ■ f*“  nr‘-n1**-i  ‘hip  tirnTrlinf  fapreetwtaliiiii 

of  the  way  in  which  manoeuvres  are  | Trims  < * F 

An  attempt  is  made  to  construct  a mmM  afhe  piles*!  actions  luViMtu.  the  i 
of  the  situation,  the  forecmt  of  the  motion  of  Se  halC^tarX^Smld^1 
the  flight  parameters. 
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The  digital  simulation  technique  is  used  flg  ensuring  greater  freedom  of  choice  of  the  com- 
putation model. 

By  dividing  the  period  of  the  manoeuvre  into  smell  time  intervals,  the  solution  of  the  equation 
of  motion  of  the  helicopter  Is  corrected  for  eMh  successive  step,  so  that  a set  of  data  describing 
the  manoeuvre  (usually  for  the  end  of  the  flight  traJecMry)  is  obtained  without  affecting  the 
manoeuvre  constraints. 

As  the  computation  technique  is  improved  the  usefulness  of  complete  simulation  of  real 
physical  processes  becomes  greater,  thus  enabling  more  detailed  analysis,  reduction  of  expensive 
experimental  research  and,  in  the  case  of  aircraft  applications,  analysis  of  the  limiting  cases, 
experimental  investigation  of  which  is  dangerous. 

Such  a simulation  is  partictdarly  justified  for  heboapters,  the  description  of  the  flight  of 
a helicopter  requiring  many  data,  but  may  also  be  of  nee  for  the  dynamic  behaviour  of  other 
controlled  objects  (such  as  airplanes,  gliders  and  mdMvcars). 
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